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Integrated photonics has enabled much progress towards quantum technologies. Many ap¬ 
plications, including quantum communication, sensing, and distributed and cloud quantum 
computing, will require coherent photonic interconnection between separate chip-based sub¬ 
systems. Large-scale quantum computing systems and architectures may ultimately require 
quantum interconnects to enable scaling beyond the limits of a single wafer and towards 
"multi-chip” systems. However, coherently interconnecting separate chips is challenging due 
to the fragility of these quantum states and the demanding challenges of transmitting pho¬ 
tons in at least two media within a single coherent system. Distribution and manipulation 
of qubit entanglement between multiple devices is one of the most stringent requirements of 
the interconnected system. Here, we report a quantum photonic interconnect demonstrating 
high-fidelity entanglement distribution and manipulation between two separate chips, im¬ 
plemented using state-of-the-art silicon photonics. Path-entangled states are generated and 
manipulated on-chip, and distributed between the chips by interconverting between path¬ 
encoding and polarisation-encoding. We use integrated state analysers to conhrm a Bell-type 
violation of 5'=2.638±0.039 between two chips. With improvements in loss, this quantum 
interconnect will provide new levels of flexible systems and architectures for quantum tech¬ 
nologies. 


Further progress towards quantum communicatiorP^^, 
sensing and computin^^wih greatly beneht from a "quan¬ 
tum photonic interconnect": an inter-/intra-chip link— 
e.g. in optical-hbre or free-space—capable of coherently 
distributing quantum information and entanglement be¬ 
tween on-chip sub-systems within a single complete quan¬ 
tum system. The signihcance of quantum interconnect 
was hrst highlighted by Kimhl^, and here we study a 
chip-based interconnect solution, which will be essential in 
many future applications and provide substantial architec¬ 
tural flexibility. Secure quantum key distribution and quan¬ 
tum communicatiorPtH, and distributed ad cloud quantum 
computin^^H^^ require interconnected on-chip subsystems 
for practical implementations. Precise quantum sensing will 
gain flexibility and versatility from on-chip generation and 
detection of entanglement, with the interaction with sam¬ 
ple performed in a different media or location, e.g., chip, 
optical hber and free spacd^^^^. Quantum computing will 
beneht from quantum interconnects through architectural 
simplihcations^^Hi^, easier i ntegr ation of materials and plat- 
forms optimised for sourc^I®^, circuhP^MUl^ qetectoi®IlII 
and many otherJ^^^ performance; and the inclusion of off- 
chip optical delays or memories. Ultimately, large-scale in¬ 
tegrated quantum systems may even exceed the area of a 
single wafer or require interconnects for architectural rea¬ 
sons. 

A quantum photonic interconnect must maintain coher¬ 
ent transmission of the qubit state a |0) + /3|I) between 
subsystems; a signihcant difference compared to the classi¬ 


cal optical interconnect in which transmitted state are either 
0 or I, and in which the relative phase is not maintained^. 
The quantum interconnect must also be capable of coher¬ 
ently interconverting between the preferred enco dings in the 
platforms and media through which it connects®^. Per¬ 
haps, the most demanding requirement for an interconnect 
is the preservation of high-hdelity entanglement through¬ 
out any manipulation, inter conversion and transmission pro¬ 
cesses within the full system. Distributing entanglemenfP^ 
between integrated chips is a key requirement and a major 
technical challenge due to the highly fragile nature of entan¬ 
glement and the potential for decoherence of quantum states 
transmitted between different chips. Path-encodin^^^ll^ —a 
photon across two waveguides—is a natural choice of ro¬ 
bust on-chip encoding for quantum in ferniat ion processing, 
however, polar isatiorPIISl^ spatial-mod^^^^^, or time-birP^ 
encoding is typically more suitable in hbre and free space 
for quantum information transmission and distribution. Al¬ 
ready there have been demonstrations of important fea¬ 
tures of quantum interconnect components, inclu ding on- 
chip entangleme nt ge neration and manipulatiorP^lf^®^^^, 
photon detectioiP^l^, interfacing of lig hUs d ifferent degrees 
of freedonJ^^^ni, and multi-chip linkJ^^^^. However, to 
date there has been no demonstration of a complete quan¬ 
tum photonic interconnect system capable of coherently dis¬ 
tributing and manipulating qubit entanglement across two 
or more integrated quantum circuits. 

Here, we demonstrate a high-hdelity quantum photonic 
interconnect. Telecom-band entangled photons are gener- 
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Figure 1: Quantum photonic interconnect and entanglement distribution between two integrated Si photonic chips, a, 

Chip-A comprises path-entangled states generation, arbitrary projective measurement A{0az,0ay), and path/polarisation conversion 
(PPG), b, Chip-B includes a projective measurement B{0bz^0by) and PPG. On the chip-A, signal-idler photon-pairs are created 
in the spiralled waveguide single-photon source. Bell states are produced when Oss is controlled to be 7r/2 or tt. Idler qubit 

initially encoded in path are coherently coupled to polarisation-encoding and transmitted through a 10m single-mode optical fibre 
(SMF), and reversely converted back to path-encoding on the chip-B. Signal qubit is analysed using A(0az, Oay) on chip-A and idler 
qubit is analysed using B{0bz^0by) on chip-B. The 2D grating coupler, behaving as the path/polarisation converter (PPG), is used 
to coherently interconvert photonic qubit between polarisation and path-encoding. Optical microscopy images of, c, the photon-pair 
source, d, the arbitrary state analyser (Inset shows the MMI splitter), and e, the 2D grating coupler PPG structure. 


ated, manipulated and distributed between two integrated 
silicon photonic chips linked by an optical-fibre. These chips 
were fabricated using state-of-the-art technologies from sil¬ 
icon photonics to enable and monolithically integrate all 
of the capabilities required for the quantum interconnect. 
Path-entangled Bell states are generated and manipulated 
on-chip. These entangled-states are distributed across two 
chips, by transmitting one qubit from one chip to the other 
via the fibre. Coherence is preserved by inter converting 
between path-encoding on chips and polarisation-encoding 
in fibre using a two-dimensional grating coupleJ^. Each 
qubit is analysed in the respective chips using thermal phase 
shifters to form integrated state analysers. We show a high- 
fidelity interconversion from path-encoding on one chip to 
polarisation-encoding within the fibre, and back to path¬ 
encoding on the second chip. We implement a rigorous test 
of entanglement—confirming a strong Bell-type inequality 
violation of 16.4cr or 15.3cr—and demonstrate the quantum 
interconnect. Together with further improvement in loss, 
this approach will provide new quantum technologies and 
applications that rely on or benefit from quantum photonic 
interconnects. 

Figure 1 shows the schematics of a chip-to-chip quantum 
photonic interconnect, generating path-entangled states on 
chip-A and coherently distributing one entangled qubit to 
chip-B, via a 10 m single mode optical fibre link. Chip- 
A and chip-B respectively have an effective footprint of 
1.2X0.5 mm^ and 0.3x0.05 mm^. On chip-A, a signal 
(As^l550.7 nm) and idler (Ai^l560.3 nm) photon pair is 
generated via the elastic scattering of two photons from a 
bright continuous-wave pump field (A^^ 1555.5 nm) inside 
2cm spiralled waveguide sources (Fig.lc), by using the spon¬ 


taneous four-wave-mixing (SFWM) nonlinear effecfP^. The 
pump is split across two sources using a multimode inter¬ 
ference (MMI) beam splitter with a near 50/50 splitting 
ratiJ^ (Fig.Id). The photon pairs are produced in either 
the top or bottom waveguides, yielding a photon-number 
entangled stat^as (11*1*) |0*0j) - |0*0j) \lsli))/'/2, 

where 6^55 is a thermally-controlled phase after the sources. 
These photons are probabilistically separated by two de¬ 
multiplexing MMIs and post-selected by two off-chip fil¬ 
ters. producing the maximally path-entangled Bell states 
1^) = (|0)s |0)i i |l)s with a 25% success proba¬ 

bility, when Oss equals to {n + l/ 2 ) 7 r or nir for an integer 
n. Subscript s and i represent the logical states of signal 
and idler qubits (more details see Supplementary Informa¬ 
tion). Then, we use an on-chip path/polarisation converter 
(PPG) to coherently interconvert the idler qubit between its 
path and polarisation-encoding. On chip-A, path-encoded 
qubit is converted to polarisation-encoded before transmit¬ 
ting across the fibre. Chip-B reverses this process, convert¬ 
ing the polarisation-encoded qubit back to on-chip path- 
encoded qubit, by using a PPG. This PPG enables entangle¬ 
ment preservation throughout the chip and fibre platforms. 
Signal and idler qubits are manipulated and measured in¬ 
dependently on two chips using arbitrary single qubit mea¬ 
surement stages A{0az^ ^ay) and B{0bz^^by)i which phys¬ 
ically consists of an on-chip Mach-Zehnder interferometer 
with an additional thermal phase shifter (Fig.Id). 

We first discuss the coherent interconversion of path and 
polarisation-encoding by using the PPG. In silicon quan¬ 
tum photonics, transverse-electric (TE) mode is usually in 
use, owing to its stronger wave guid e confinement and conse¬ 
quently enhanced SEWM effecfP^. Thus, produced photons 
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must be guided iu TE-modes, which is easily achieved by iu- 
jectiug pump light iu this mode usiug a ID TE-gratiug cou¬ 
pler. Our PPG is implemeuted usiug a 2D gratiug coupler 
(Eig.le), where TE-polarised light comiug from two uearly 
orthogoual waveguides is c ombiu ed iuto two orthogoual po¬ 
larised compoueuts of lighlP^^^. Iu this way, the polarisa- 
tiou states of photous received by the fibre is determiued by 
the two-waveguide ou-chip states, aud vice versa. This pro¬ 
vides a cohereut iutercouversiou betweeu path-eucodiug aud 
polarisatiou-eucodiug. Details are provided iu the Methods 
aud SI. To coufirm the PPG cohereut mappiug, we prepared 
arbitrary bright-light polarisatiou states usiug bulk opti¬ 
cal compoueuts aud coupled them iuto the ou-chip receiver 
(Eig.2a). The 2D gratiug coupler couverted the polarisa¬ 
tiou states iuto path-eucoded states, which were theu aual- 
ysed ou-chip by implemeutiug a full state tomograph}®. 
We prepared a set of six polarisatiou states ppoi^ aud mea¬ 
sured the correspoudiug ou-chip path states ppath] these 
states are showu as Bloch (or Poiucare) vectors iu Eigures 
2b aud 2c, respectively. The distauce betweeu the states 
cau be described by the state fidelity, which is defiued as 
Estate = ■ ppath ■ ^/^])^■ The mean fidelity of 

the six measured states is 98.82±0.73%. Au example of a 
recoustructed deusity matrix for the path-eucoded state |+) 
is showu iu Eig. 2d (full data are provided iu Eig.S5). Theu, 
we fully quautified the PPG process usiug a quautum pro¬ 
cess tomograph}®. This is mathematically described by a 
process matrix y, defined by ppath = T.mni^niPpoiElxmn), 
where Ei are the Ideutity matrix I aud Pauli matrices 
X, T, aud Z, respectively. By subjectiug the ppoi states 
iuto the PPG aud measuriug the ppath states, we deter¬ 
miued the process matrix y of the PPG, showu iu Eig.2e. 
We hud a high process fidelity of 98.24±0.82%, defiued as 
^process = Tr[xideai • x], where Xideai is the ideal process 
matrix with Xideai[I^I] = 1 - ^5 aud Z amplitudes of the 

measured matrix y represeut the probability of a bit-flip or 
phase-flip error ou the PPG iutercouversiou. 

A ^50 mW GW pump was iujected iuto ou-chip sources 
to create photou pairs. Photous were detected usiug two su- 
percouductiug uauowire siugle-photou detectors (SNSPDs) 
with ^50% efhcieucies aud ^800 Hz dark couuts. Goiuci- 
deuces were recorded usiug a time iuterval aualyser. After 
the chip-A, a meau rate of 500—800 Hz photou pairs was 
observed, while after the two chips we obtaiued 8—12 Hz 
meau coiucideuces. The pump light propagates colliuearly 
with siugle photous, aud this allows a closed feedback loop 
to track photous aud mouitor state stability throughout the 
chip-to-chip quautum iutercouuect system. Details are pro¬ 
vided iu the Methods aud SI. 

We uext coufigured chip-A to produce eutaugled states. 
Sigual aud idler photous were respectively collected at ports 
D1 aud D2 of the chip-A, aud routed to SNSPDs. Goutiu- 
ually scauuiug Oss^ we observed the A-classical iuterfereuce 
aud A/2-quautum iuterfereuce friuges with a visibility {V = 
1 — Njain/Nmax) of 99.99±0.01% aud 99.36±0.17%, respec¬ 
tively (Eig. 3a). The high visibility of this double-frequeucy 
friuge is a siguature of high-quality photo u-uumber eutau- 
glemeut produced iuside the chip-A^^^l^. The high vis¬ 
ibilities arise from well-balauced MMI splitterJ^ aud a 
good spectral overlap betweeu two photou-pair source^^^^. 
Theu, the photou-uumber eutaugled state evolves iuto the 
path-eutaugled Bell states, |4>)+ or |4>)“, by settiug Oss to 


a a|H>+j8|V> 


a|0>+j8|1> 



b "Polarisation Bloch" 

\H) 


C "Path Bloch" 

|o> 



d state matrix p of the | +> 



6 Process matrixy of the PPC 



Figure 2: Interconversion between polarisation-encoding 
and path-encoding, a, Initial arbitrary polarisation-encoded 
states a\H) -h/3|E) {\H) and \ V) are two orthogonally polarised- 
states) were prepared by using a set of polariser (P), half-wave 
plate (HWP) and quarter-wave plate (QWP). A fibre polari¬ 
sation controller was used to compensate polarisation rotation 
in the fibre. The PPG converted polarisation states into path- 
encoded states o|0) -h /3|1), where |0) and |1) denote path states 
in two waveguides. The path-encoded states were analysed us¬ 
ing B{0bz,0by) to implement state tomography, b and c, the 
Bloch sphere representation of ideal polarisation-encoded states 
{\H), \ V), \D), \A), \R), \L)} in bulk optics, and measured path- 
encoded states {|0), |1), |-h), |—), I + ^), I — i)} on chip. Indicated 
fidelity represents the mean over all six states, d, Reconstructed 
density matrix of the |-|-) path-encoded state corresponding to 
the \D) polarisation-encoded state, e, Reconstructed process 
matrix y of the PPG using the quantum process tomography. 


7r/2 or TT. The entangled-qubits were separated and coher¬ 
ently distributed across chip-A and chip-B using the PPG 
iutercouversiou. We measured correlation fringes across the 
two chips, by collecting signal photons at port D1 on chip-A 
and idler photons at port D3 on chip-B, and simultaneously 
operating A{0azi ^ay) and B{0bz^ ^by) on two chips. Eig¬ 
ures 3b and 3c show the entanglement correlation fringes 
for the Bell states |4>)+ and |4>)“ as a rotation of Oby on 
chip-B, with Oay on chip-A set at {0, 7r/2 , tt, 37 r/2}. These 
results are in good agreement with their theoretical predic¬ 
tions of cos^\{0 ay ~ and cos^\{0 ay T 

The fringes exhibit a mean visibility of 97.63±0.39% and 
96.85±0.51%, respectively, which is far beyond the critical 
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visibility required of 1/v^ to violate the Bell inequalit}®. 
We now have shown that a very high-quality entanglement 
is produced on the chip-A, and distributed over the fibre to 
the chip-B, coherently interconnecting the two chips. 

To more strictly verify the existence of entanglement 
across the two chips, we directly measured the Bell-CHSH 
(Clauser-Horne-Shimony-Holt) inequalit}?^, defined as: 


5 = ||(Ai,5i) + (Ai,52) + (A2,5i) - (A2, 52)11 < 2 (1) 

where Ai and Bi briefly denote the projectors A(0, Oay) and 
5(0, Oby) on the two chips. Correlation coefficients (A^, 5^) 
were measured, when Oay on chip-A was set to {0,7r/2} and 
Oby on chip-B was set to {7r/4, 37r/4}. Full data of (A^, Bi) 
is provided in Fig.SG. Using equation (1), we obtained the 
directly measured Schsh parameters of 2.638A0.039 and 
2.628±0.041 for the two Bell states |4>)+ and |4>)“, respec¬ 
tively. These Schsh violate the Bell-CHSH inequality by 
16.4 and 15.3 standard deviations, strongly confirming that 
the two photons after distributed across chip-A and chip-B 
are highly entangled, and therefore verifying a high-quality 
quantum photonic interconnect between two chips. In addi¬ 
tion, we estimate the maximally achievable Sfringe parame¬ 
ters of 2.761T0.011 and 2.739T0.015 for the |^)+ and |^)- 
states, from the mean visibility of the correlation fringes 
(shown in Fig.3) according to Sfringe = 2\/2VE^. Figure 
4 illustrates a good agreement of the Schsh and Sfringe 
parameters. 

We demonstrate high-fidelity entanglement generation, 
manipulation, inter conversion, distribution and measure¬ 
ment across two separate integrated photonic circuits, suc¬ 
cessfully demonstrating the first chip-to-chip quantum pho¬ 
tonic interconnect. Path-polarization interconversion pre¬ 
serves quantum coherence across the full interconnected 
chip-fiber-chip system, demonstrating beyond single-chip 
implementation of a quantum photonic experiment. The 
efficiency of this interconversion process can be further 
improved by engineering the geometry of grating cou¬ 
pler or utilising other polarisation control technique d^^l^^ l 
Other inte rconve rsion approaches, e.g., path-orbital angular 
momenturrP^^^ or path-time binJ^ may further enrich this 
quantum photonic interconnectivity. The use of silicon— 
allowing large-scale integratioiP^and compatibility wi th mi- 
croelectronics and telecommunications infrastructurd^^^, 
benefiting fr om th e classical optical interconnect technol¬ 
ogy on silicorP^^, and also o ffering abili ty to monoli thical ly 
integrate photon sourcd^^^, circuifP^^^ and detectoJ ^^ * ^^ * — 
would position this quantum photonic interconnect technol¬ 
ogy at the heart of building practical, robust and scalable 
silicon-based quantum hardware for future and networks. 
This work opens the door to multi-chip integrated quantum 
photonic technologies, which would be capable of processing 
quantum information on extremely small and stable chips 
and also capable of robustly distributing and transmitting 
quantum information between chips. 


Methods 

Devices design and fabrication. The devices were fabri¬ 
cated on the standard silicon-on-insulator wafer with a 220 nm 
silicon layer and a 2 buried silica oxide layer. MMI couplers 
were designed as 2.8 pmx27 pm to get a balanced splitting ratio 


a Vj, = 99.99±0m% V\/2 = 99.36±0.17% 



055 


b »/[|0>q = 97.63±O.39% 



C »/[|0>-] = 96.85±O.51% 



Figure 3: Entanglement fringes, a, A-classical interference 
(cyan) and A/2-quantum interference (red) fringes measured on 
chip-A. Bright-light was measured (normalised) at port Dl, and 
coincidences were collected (accumulated 20s) between ports Dl 
and D2. Oss was rotated to produce fringes, as A{0az,0ay) was 
set as a Hadamard gate. Photons are bunched or anti-bunched 
when Oss is utv or (n -h l/2)7r. b and c. Entanglement correla¬ 
tion fringes for the Bell states |4>)^ and |4>)“ after distributed 
across the two chips. Coincidences were collected (accumulated 
30s) between ports Dl and D3. Oby on chip-B was continually 
rotated {Obz = 0) to obtain the fringes, as A{0az, Oay) on chip- 
A was projected onto {|1), |0), |—), |-h)} basis by setting Oay to 
{0,7r/2, TT, 37r/2} and Oaz to 0. The indicated visibility repre¬ 
sents the mean over all four fringes. Error bars are given by 
Poissonian statistics. Accidental coincidences are subtracted. 


(Eig.ld). MMIs offer a large bandwidth and a large fabrication 
tolerance. Spiralled waveguide sources with a 2-cm length were 
used to create photon-pairs. The ID grating couplers consist of 
a periodic 315 nm silicon layer with a 630 nm pitch. The 2D 
grating couplers include 10 pmxlO pm hole arrays with a 390 
nm diameter and a 605 nm pitch. Resistive heaters with a 50 
pm-length were designed and formed by a Ti/TiN metal layer. 
The devices were fabricated using the deep-UV (193 nm) lithog¬ 
raphy at LETI-ePIXfab. Silicon waveguides were 220 nm fully 
etched, while ID and 2D grating couplers were 70 nm shallow 
etched. The devices were covered by a 1.6 pm silica oxide layer. 

Devices characterisations. Optical accesses and electric 
accesses were independently controlled on two chips (Eig.Sl). 
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Figure 4: Verification of chip-to-chip entanglement distri¬ 
bution and quantum photonic interconnect. The two Bell 
states were distributed across two Si chips. The S parameters 
were obtained using two approaches. Green dotted columns are 
the directly measured Schsh using Eq. (1). Pink columns are 
the maximal achievable Sfringe, estimated from the mean visi¬ 
bility of correlation fringes in Fig.3. The Schsh and Sfringe pa¬ 
rameters are in good agreement. Black and blue dashed lines de¬ 
note the classical and quantum boundary. These results confirm 
the high levels of entanglement after distributed across chips, 
and the high fidelity of quantum photonic interconnect. Co¬ 
incidences of each measurement were accumulated for 60s and 
accidental coincidences are subtracted. Error bars (±1 s.d.) are 
given by Poissonian statistics. 


Optical accesses were achieved using V-groove single modes fibre 
arrays with a 127 pm pitch. Fibres were titled with an angle of 
10~12 degrees to guarantee both grating couplers work at the 
required wavelengths. Excess loss of ID and 2D grating couplers 
were about -4.8 dB and -7.6 dB at peak wavelengths, respec¬ 
tively (Fig.S2). Extinction ratio of ID and 2D grating couplers 
were measured to be larger than 20 dB and 18 dB, respectively. 
We estimated losses from different contributors in the system: 
-6 dB from off-chip filters, -6 dB from SSNPDs, -9~9.5 dB from 
ID grating couplers, -15~15.5 dB from 2D grating couplers, -6 
dB from demultiplexing MMIs, and -8^9 dB from MMIs excess 
loss and propagation loss in waveguides. Totally, signal and idler 
photons respectively experienced -36~38 dB and -18~19 dB at¬ 
tention. The fibre channel interconnected two chips has -15~16 
dB attention mainly from the 2D grating couplers, which can 
be further improved by engineering the geometry of coupler and 
waveguide. 


All thermal-driven phase shifters were controlled using home¬ 
made electric controllers. Wire bounding technology was used to 
contact heaters’ transmission lines. Optical power was recorded 
as a function of electric power added on heaters. The optical— 
electric power contour was fitted and used to construct the map¬ 
ping between the required states and electric power. Fig.S3 
shows calibration results of chip-A’s and chip-B’s state analy¬ 
sers. To avoid the influence of temperature variation, both chips 
were mounted on temperature stabilised stages. Fibre alignment 
was automatically recoupled using piezo-electronic stacks. Fig.S4 
shows the stability of the chip-to-chip system, which were main¬ 
tained constant more than 30 mins. This indicates path-encoded 
states on the two chips are very stable and polarisation-encoded 
states in the fibre channel are also well-stabilised. 
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Appendix A: Quantum state evolution 

A pair of signal (Ag) and idler (A^) photons are created 
in the spiralled waveguide source by annihilating two pump 
photons (Ap), based on the x(3) SFWM nonlinear optical 
effect [1, 2], whose Hamiltonian is approximately described 
as: 

Hocalapa] + agO^^ai (Al) 

where a|, aj, and a\ are the creation operators, and a^, 
Up, and ai are the annihilation operators, for the involved 
signal, pump, and idler photons. The CW bright pump-light 
is equally split into two spiralled sources to produce a pair 
of signal and idler photons, which are coherently bunched 
either at the top or bottom waveguides as [3, 4]: 

10) = (|l,li), |0,0,), - |0,0,), |l.li),)/V2 (A2) 

where 0 and 1 represent photon number at the top and 
bottom waveguides with subscript t and 6, respectively. Oi 
and Os are phases of the signal and idler photons controlled 
by the Oss phase shifter. The difference between Oi and Og 
is negligible and the state can be simply rewritten as the 
typical NOON states format: 

|0) = (|2)jO),-e2^-MO)j2),)/V2 (A3) 

Then, signal and idler photons are split by the demulti¬ 
plexing MMI couplers. Above photon-number entanglement 
state evolves to: 


1 ^) ~ 

—i 

+ 

+ 

—i 

+ | 0 s 0 i)tt | 0 s 0 i)t 6 

+ |0s0i)tt |0s0i)t6 bb]} 

And after the crossing, the state can be rewritten as: 


1 ^) ~ | 0 s 02 )t 6 

—i lOgO^)^^ 

+ |0sli)tt |0s0i)t6 \^s^'i)bt 
+ \^s^i)u |0s0i)t6 \^s^i)bt |0s0«)66] 

—i |0s0i)^^ |0s0i)^^ 

+ | 0 s 0 i)tt | 0 sli)t 6 

+ |0s0i)tt |ls0i)t6 \^s'^i)bb]} 


The signal and idler photons are respectively selected by 
two off-chip dense-wavelength-demultiplexer (DWDMs) fil¬ 
ters, and following items of the state are post-selected with 
a 25% probability: 

10 ) = ~^[\'^s^i)tt\^s^i)tb\^s'^i)bt\^s^'i)bb 

|0A)« |lA)*jO,Oi),jO,li),J 

(A4) 

Reform the state using the logical qubit representation 
as: 

l0) = [|O),|O),-e2*^-|l)Jl)J/x/2 (A5) 

where |0)s (|0)i) and |l)s (|l)i) denote signal (idler) pho¬ 
ton’s path states in two waveguides. When Oss is chosen to 
be (n+l/2)Ti or mz for an integer n, we respectively obtain 
the two Bell states: 

|$)^ = [|0)JO),±|l)Jl)J/y2 (A6) 


Appendix B: Devices characterisation 

Figure SI shows the experimental setup for the distribu¬ 
tion of entanglement between two integrated silicon pho¬ 
tonic devices. Bright light around a wavelength of 1555.5 
nm was collected from a tunable laser (Tunics-BT) and fur¬ 
ther amplified using a high-power EDFA (Pritel). The am¬ 
plified spontaneous emission (ASE) noise was suppressed 
using a DWDMO (Openti) filter. Next, the pump light was 
injected into the chip-A by using an 8-channel single-mode 
fibre array (OZ-Optics) with a 127 pm pitch and a 10 de¬ 
gree polished angle. A fibre polarisation controller (PC) was 
used to guarantee that TE-polarised light was launched into 
the device through the ID TE grating coupler. Entangled 
photon-pairs were produced on the chip-A, and coherently 
distributed to the chip-B via a 10-m single-mode fibre. Po¬ 
larisation rotation in the fibre channel was compensated by 
using another PC. Two off-chip DWDMs with a 200 GHz 
channel space and 1 nm IdB-bandwidth were used to sep¬ 
arate signal and idler photons. We selected the photons 
which are equally 3-channels away from the pump, that is 
Ap — As =Xi — Xp = 4.8nm. This gave a high extinction ra¬ 
tio to efficiently remove the remained pump from photons. 
Moreover, DWDMl after the chip-A only picked up the sig¬ 
nal photons with As= 1550.7 nm, while DWDM2 after the 
chip-B only selected the idler photons with A^= 1560.3 nm. 
Note DWDMs induced ~3dB loss for each photon. Pho¬ 
tons were detected using two fibre-coupled superconducting 
nanowire single-photon detectors (SNSPDs) mounted in a 
closed cycle refrigerator. Polarisation before SNSPDs was 
optimized to maximize detection efficiency up to ~50%. Co¬ 
incidences were recorded by using a time-interval-analyser 
(PicoHarp 300) in a 450 ps integration window, which was 
chosen according to the jitter-time of SNSPDs. 
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EDFA: Erbium-doped fibre amplifier 
DWDM; dense wavelength demultiplexer 
PC: polarisation controller 
PEC; piezo-electronic controller 
PD; photodetector 

SNSPD: superconducting nanowire single-photon detector 
TIA; time-interval-analyser 



Figure 5: Experimental setup for the distribution of entanglement between integrated Si photonic chips, chip-A and chip-B. 


1. Grating couplers characterisation 

Figure S2 shows the measured spectrums of the ID and 
2D grating couplers. Peak wavelengths of both gating cou¬ 
plers are dependent on the angle between fibre array and 
chip, and they are both around 1555.5 nm when the rel¬ 
ative angle is in the range of 10—12 degrees. Excess loss 
of ID and 2D grating couplers is about -4.8 dB and -7.6 
dB at the peak wavelengths, and their IdB-bandwidths are 
around 27 nm and 30 nm, respectively. The ID grating cou¬ 
plers consist of a periodic 315 nm silicon layer with a 630 
nm pitch. The 2D grating couplers include 10 pm x 10 pm 
hole arrays with a 390 nm diameter and a 605 nm pitch 
(Insets of Fig.S2). Optimised angles for the two chips are 
slightly different, owing to the wavelength difference of sig¬ 
nal and idler photons and also fabrication deviation of the 
devices. Loss of grating couplers can be further reduced by 
engineering the grating structure and positioning reflection 
mirrors under the grating [5, 6]. Note that other on-chip 
polarisation control and diversity devices can be explored 
for quantum linking between chips using entanglement [7]. 


2. Tomography stages characterisation 

The two chips were fixed on two copper PCBs and 
thermal-heaters were wired-bounded to electric pads on the 
PCBs. Home-made computer-interfaced heater drivers were 
used to independently control all heaters on the two chips 
(Figure SI). The output optical power was recorded as a 
function of electric power added on heaters, from which the 
relationship between states and electric power was recon¬ 
structed. A least-squares minimisation algorithm was used 
to fit the 0-E (optical power and electric power) contour 
and find responding powers for different states. Figure S3 
shows the calibration results of chip-A’s and chip-B’s projec¬ 
tors, A{0az,0ay) and B(0bz,0by)^ by simultaneously scan¬ 
ning Oay and Oaz or Oby and Obz phase shifters, respec¬ 
tively. To calibrate chip-B’s B{0bz,0by), the input state 
needs to be known in advance otherwise it is difficult to 
access phase offset of the Obz phase shifter. We used the 
ID TE-grating coupler as an on-chip polariser to guaran¬ 
tee TE-polarised state was injected. Then we kept the fi¬ 
bre untouched and smoothly switched the input state to 
the 2D grating coupler on the same chip. This state is an 



Wavelength (mn) 

Eigure 6: Spectrums of the ID and 2D grating couplers. Peak 
wavelength with maximal transmission is determined by the ti¬ 
tled angle of fibre. These spectrums were measured when the 
relative angle was optimised to be 10-12 degrees. In this case, 
both grating couplers work near 1555.5 nm peak wavelengths. 
Insets show the optical microscopy images of the grating cou¬ 
plers. 


anti-diagonal state for the 2D grating coupler. Then we 
determined all phase information of chip-B’s B{0bz,0by)- 
Similarly we calibrated chip-A’s A{0az,0ay)- 


3. Systematic stability measurement 

A classical reference frame is necessary for real-life scenar¬ 
ios. The phase matching condition of the SFWM processing 
allows the generated photons propagate collinearly with the 
pump light. Then, we can use the pump light to close a 
feedback loop to track single photons on both two chips 
and also in the optical fibre channel, and to keep state sta¬ 
bility in the chip-to-chip system. The demultiplexing MMIs 
on chip-A split the pump into two parts, chip-A’s projec¬ 
tor and the chip-B. Half was used to feed-forwardly track 
photons at the chip-A side; the other half was transmitted 
through the fibre channel together with the idler photons, 
and further used to track photons in the optical fibre and 
the chip-B side. Figure S4 shows the stability of the chip- 
to-chip system. Fibre alignment was feed-forwardly recou¬ 
pled each 1 min by using a piezo electronic stacks (PEC). It 
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a Chip-A Data b Chip-A Fit C Chip-B Data d Chip-B Fit 



0AY power (mVA) 0 ay power (mVA) 0 by power (mVA) 0 by power (mVA) 


Figure 7: Calibrations of chip-A’s and chip-B’s projectors. Optical power was recorded as a function of electric power added on 
heaters, (a) and (b) are experimental data and fit of A{0az,0ay) and (c) and (d) are experimental data and fit of B{0bz^0by) in 
which axes are electric power added on the Oay and Oaz or 6by and Obz heaters, and contour color represents normalised optical 
power. This optical-electric contour was fitted by using the least-squares minimisation algorithm to find the demanded states. White 
dots at crossers denote the searched projective states {|0), |1), |+), |—), | + ^), | — i)} in the map of electric power. 



Figure 8: Stability of the chip-to-chip system. Curves represent 
the measured optical power and self-normalised, as a function of 
time. Purple one is obtained at chip-A’s port Dl, and pink one is 
measured at chip-B’s port D3. A(0az, Oay) and B{0bz, Oby) on 
both two chips were set as the |+) projectors. Fibre alignment 
was feed-forwardly recoupled each Imin. 


the state fidelities defined as the overlapping between ppoi 
and ppath- By subjecting polarisation-encoded states into 
the PPG and measuring output path-encoded states (deter¬ 
mined using state tomography), we determinate the process 
matrix x of tho PPG processing for a fixed set of opera¬ 
tors Ei^ where Ei, E 2 , and Es are simply chosen to be the 
Pauli operators and Eq to be the Identity operator. Twelve 
parameters in the PPG’s process matrix x need to be de¬ 
termined. Four ideal polarisation states ppoi, for example, 
{|F7), |i/), \D), 1^)}: were chosen as the input, and recon¬ 
structed path-encoded states ppath {lO), ll), 1+), I + 
were used as the output. The process matrix x is described 
explicitly as [8]: 

X = AQA (Cl) 

where the matrix Q is determined by the reconstructed den¬ 
sity matrix ppath of |0), |1), |+), | i) and the matrix A is 
described as: 


shows that path-encoded states on the two chips are stable 
and polarisation-encoded states in the fibre channel are also 
well-stabilised. A longer-time stability maintenance would 
require an active polarisation control before the chip-A chip 
and also between two chips. 


X = 


I X 
X -I 


(C2) 


Appendix D: Entanglement correlation coefficients 


Appendix C: State and process tomography 

Initial polarisation encoded states, {\H), \V)^ |T)), |A), 
|i7), 11/)}, were prepared in bulk optics and injected to a 
copy of the chip-B. The converted path-encoded states, {|0), 
|1), |+), |—), I + i), I — were analysed on-chip by im¬ 
plementing state tomography [8]. Each state was measured 
in six basis states {|0), |1), |+), |—), | + i, | — '<^)}- A max¬ 
imum likelihood search algorithm was used to reconstruct 
the most likely legitimate state from measurement results. 
Figure S5 shows the ideal polarisation-encoded states ppoi 
and the six reconstructed path-encoded states ppath ^ and 


Gorrelation coefficients (Ai^Bi) or {Ai{0AY)^ Bi{0BY)) 
inside the Bell-GHSH inequality are defined as the nor¬ 
malised correlation value of the measured coincidences data, 
shown as below [9, 10], equation (Dl). In our experiment, 
we respectively measured the coincidences with a chip-A’s 
Oay set of {0, ti/2, ti, 3ti/2} and a chip-B’s Oby set of 
{ti/4, 3ti/4, 5ti/4, 7ti/4}, to get the correlation coefficients 
of (A(0),B(7r/4)), (A(0),B(7r3/4)), {A{tt/2), B{tt/4)), and 
(A(7r/2), 5(7r3/4)). Figure S6 shows the measured correla¬ 
tion coefficients of the two Bell states. These correlation 
coefficients were consequently used to calculate the Bell- 
GHSH S parameter using equation (D2). 
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Figure 9: Density matrix of ideal polarisation-encoded states and reconstructed path-encoded states, (a)-(f). Ideal density matrix 
Ppoi of the polarisation-encoded states, |i7), |y), |D), |A), \R)^ and|L), prepared in bulky optics, (g)-(i). Reconstructed density 
matrix ppath of the six path-encoded states, |0), |1), |H-), |—), | + i) and | — ^), measured on chip. The state fidelities between ppath 
and ppoi are estimated to be 98.66%, 98.38%, 99.34%, 99.44%, 97.66%, and 99.46%, respectively. 



Figure 10: Measured correction coefficients of two Bell states |4>)^ and |4>)“ after distributed across the chip-A and chip-B chips. 
Coincidences of each measurement were accumulated for 60s. Accidental coincidences are subtracted for all data. Standard deviation 
of the correction coefficients are calculated from an evolution of Poissonian photon statistics. 


IA(0 'I B(0 )] — +tt,0by + tt) - C{9ay,(^by + tt) - C{0ay + tt, ^gy)] , , 

[C{0ay,Sby) + C{9ay + Tr,0BY + tt) + C{9ay,()by + tt) + C{9ay + tTjOby)] 

S =11 {Ai{9ay), Bi{9by)) + {Ai{9ay), B2{9by)) + {A2{9ay), Bi{9by)) — {A2{9 ay), B2{9 by)) || (D2) 
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